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ON THE PROBABILISTIC APPROACH TO THE SOLUTION OF 
GENERALIZED FRACTIONAL DIFFERENTIAL EQUATIONS OF 
CAPUTO AND RIEMANN-LIOUVILLE TYPE 


M. E. HERNANDEZ-HERNANDEZ, V. N. KOLOKOLTSOV 


Abstract. This paper provides a probabilistic approach to solve linear equa¬ 
tions involving Caputo and Riemann-Liouville type derivatives. Using the 
probabilistic interpretation of these operators as the generators of interrupted 
Feller processes, we obtain well-posedness results and explicit solutions (in 
terms of the transition densities of the underlying stochastic processes). The 
problems studied here include fractional linear differential equations, well an¬ 
alyzed in the literature, as well as their far reaching extensions. 


1. Introduction 

The theory of fractional differential equations is a valuable tool for modeling a 
variety of physical phenomena arising in different fields of science. Their numer¬ 
ous applications include areas such as engineering, physics, biophysics, continuum 
and statistical mechanics, finance, control processing, econophysics, probability and 
statistics, and so on (see e.g. [3], [19], (32], [36], and references therein). Fractional 
ordinary differential equations (FODE’s) and fractional partial differential equa¬ 
tions (FPDE’s) have been used as more accurate models to describe, e.g., relaxation 
phenomena, viscoelastic systems, anomalous diffusions, Levy flights and continuous 
time random walks (CTRW’s) (see, e.g., [2], [20]-[21], [25], [27], [29], [36]). 

To solve this type of equations different analytical and numerical methods have 
been investigated. Analytical methodologies include the Laplace, the Mellin and the 
Fourier transform techniques H, [13, E5-E3], and the operational calculus method 
m-m, m- Regarding the numerical approaches, we can mention the fractional 
difference method, the quadrature formula approach, the predictor-corrector ap¬ 
proach as well as some numerical approximations using the short memory principle 
amongst others (see, e.g., 0 - 12 ], 13, m and references therein). In the present 
article we focus on a probabilistic approach to solve linear differential equations 
involving natural extensions (from a probabilistic point of view) of the Caputo and 
Riemann-Liouville (RL) derivatives of order /3 G (0,1). 

The use of probability theory to solve classical differential equations is an ef¬ 
fective approach to obtain their solutions by relating them with boundary value 
problems of diffusion processes. In the fractional framework, connections between 
probability and FPDE’s have been also analyzed in the literature [H], [13-EI], E3, 
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m, [35]. For instance, the probabilistic interpretations of the Green (or fundamen¬ 
tal) solution to the time-fractional diffusion equation and the time-space fractional 
diffusion equation are already known (see references above). 

In this paper we employ similar probabilistic arguments (transforming the orig¬ 
inal problem into a Dirichlet type problem) to study linear equations involving a 
general class of Caputo and Riemann-Liouville type operators. These operators can 
be thought of as the extensions of the classical Caputo and RL fractional deriva¬ 
tives, respectively. As was shown in [^, they can be obtained as the generators 
of Markov processes interrupted on an attempt to cross a boundary point. The 
Caputo and RL derivatives of order /3 S (0,1) are particular cases arising by stop¬ 
ping and killing a ,5-stable subordinator, respectively. This fact allows one to solve 
fractional equations as particular cases of more general equations involving D— and 
—operators of the type and respectively (see definitions later). 

The problems we are addressing here are the following: 

(i) the linear equation with the Caputo type operator: 

-£'i+,w(f) = AM(f)-g(t), te{a,b], u{a) = Ua, (1) 

for a given A > 0, a bounded function and Ua G K. Since there is a rela¬ 
tionship between Caputo and RL type operators (similar to that between the 
classical Caputo and RL derivatives), we also studied the problem with the 
RL type operator: 

- Ol'^uit) = Xu{t) - g(t), tG{a,b], u(a) = 0, (2) 

(ii) the generalized mixed fractional linear equation 

d 

- ^ u{ti, ...,td) = Xu{ti, ...,td)- g{ti,.. .,td), (3) 

with some prescribed boundary condition, where denotes either the 

RL type operator — tiD^af+ or the Caputo type operator — ■ The left 

subscript ti indicates the operator is acting on the variable ti. 

Fractional linear differential equations with Caputo derivatives of order jd S (0,1) 
are particular cases of equation ©• They have been extensively investigated by 
means of the Laplace transform method. Hence, it is known that [7], [5^-[55] 

L>f+*w(t) = Xu(t) + g{t), u{a) = Ua, A S K, (4) 


for P G (0,1) and a given continuous function g on [o,6], has the unique solution 

u(t) = UaEp [ - A(t - a)'^] -f f g{r){t - rf-'^Ep^g ( - A(t - rf) dr, t G (a, b], 

J a 

( 5 ) 

where Ep and Ep^p denote the Mittag-Leffler functions (see definitions later). This 
solution can be written in terms of /3—stable densities by means of the integral 
representation of the Mittag-Leffler functions given in ()49|) . The probabilistic ap¬ 
proach introduced here will give this expression directly once one writes down the 
expectations involved in the general stochastic representation (H71) . On the other 
hand, using the results obtained here and the uniqueness of solutions, we obtain a 
pure probabilistic proof of the well-known equality in ()49p . 
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Apart from the classical Caputo derivatives, operators include, as simple 

particular cases, the multi-term fractional derivatives with 

non-negative functions uji. Hence, as another example of (ED, our approach also 
applies to the multi-term fractional equation 

k 

'^u:,{t)D^f^u{t) =-\u{t)+g{t), /3i e (0,1), t G M+, ( 6 ) 

i=0 

with some given functions g and z = 1,..., fc. The explicit solution for the case 
when LOi are constants and ([BD is a commensurate equation (i.e., the quotients Pi/ (3j 
are rational numbers for all i,j), has been analyzed by the reduction of the original 
problem to either a single- or multi-order fractional differential equation system (see, 
e g-, H, 0 , and references therein). An approximation for its solution has been 
also studied, e.g., in [7]. Our approach encompasses not only the commensurate 
case with constant coefficients uii but also the more general case with non constant 
coefficients uJi{-) and, even more generally, functions Pi(t). 

The fractional counterpart of equation ([SD is the mixed fractional equation 

- tiD^+u{ti,t2) - fPQj^^u{ti,t2) = Xu{ti,t2) - g{ti,t2), /3,aG(0,1), (7) 

subject to some boundary condition, where A > 0 and 5 is a given function on 
[0,&i] X [ 0 , 62 ]- The probabilistic approach presented here provides the explicit 
solution in terms of /3— and a—stable densities. To our knowledge this type of 
mixed fractional equations with Caputo and RL derivatives of order in (0,1) has 
not been explored explicitly in the literature. The solution to more general cases 
can be obtained by following an analogous procedure, for instance, to solve the 
equation 


_ ^ i/)Pi ^ ^ 2 ) = Xuifi, ...,td)- gih,td), ( 8 ) 

with being either the RL or the Caputo derivative. 

The probabilistic ideas behind the solution to ([1]) and (ED can be used to solve 
the linear equation with non-constant coefficients, as well as more general FPDE’s 
involving operators of the type — Ax, where — Ax denotes the generator of 

a Feller process on acting on the variable x G Some specific cases include 
the widely studied time-fractional diffusion equation (see, e.g., 0 , m, ES] and 
references above) 


d^u{x,t) d^u{x,t) 


0 < /3 < 1, C > 0, 


dtP dx"^ 

where (x, t) denotes the space-time variables. These cases will be addressed in a 
forthcoming paper in preparation. 

It is worth mentioning that equations involving and Z?—operators do not 
usually have solutions in the domain of the generators and —D^^, re¬ 

spectively. The existence of such solutions is restricted to a specific value in the 
boundary condition. Thus, as usual in classical stochastic analysis, by introducing 
the concept of a generalized solution we are able to study the well-posedness (in 
a generalized sense) for these equations. To illustrate this, consider the very-well 
known ordinary differential equation (ODE) 


u{t) = —Xu{t)-\-g{t), tG(0,6], 'u(O) = uo. 
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where 6 > 0, g S C[0,6] and A > 0, and whose solution is 

u(t)=uoe^*'+ f exp{X{t — r)}g{r)dr, t € {0,b]. (9) 

Jo 

Probabilistically, this problem can be thought of as the boundary value problem 
associated with the deterministic linear motion on (—oo, &] which is stopped at 
reaching the boundary point t = 0. In this case, the semigroup {S's}s>o of the 
deterministic process is given by Ssf{t) = fit — s) for any t S {—oo, b] and / € 
C'(—cx),&], whilst the semigroup {S'°+*}s>o for the stopped process is 5's+*/(t) = 
/(max{0,t — s}) for any t G [0,6] and / G C'[0,6]. Hence, the resolvent operator 
of the semigroup 5'°+* provides equation ([9|) as the unique solution in the domain 
of the generator (the space (^^[O, 6]) if, and only if, uq = y<?(0). Otherwise this 
solution is (in our terminology) only a generalized solution as it can be obtained 
as a limit of solutions in the domain of the generator. Moreover, in this case the 
generalized solution is also a classical (smooth) solution lying on C[0,6] O 0^(0, 6] 
instead of C'^[0, 6]. Similar situations occur when considering fractional differential 
equations. We will see that the solutions found in the literature are usually solutions 
in the generalized sense, as they usually do not belong to the domain of — or 
—as generators of Feller processes. 

The main contribution of this work relies on providing well-posedness results and 
explicit integral representations of solutions to linear equations with Caputo and 
RL type operators. We pay attention to the existence of two types of solutions: 
solutions in the domain of the generator and generalized solutions. The latter con¬ 
cept defined for rather general, even not continuous, functions g in O- Moreover, 
all solutions are given in terms of expectations of functionals of Markov processes. 
From the point of view of numerical analysis, this representation can be exploited 
to obtain numerical solutions to a variety of problems by performing Monte Carlo 
techniques. Simulation methods have been effectively used for classical differential 
equations and, in recent years, different methods for evaluating path functionals of 
Levy processes have been actively researched (see, e.g., [12]- [Hj, [23]). Finally, by 
using the monotonicity of the underlying processes, we obtain explicit solutions in 
terms of the transition densities of the Markov processes involved. 

The paper is organized as follows. The next section provides some standard 
notation and gives a quick summary about Caputo and RL derivatives as well 
as /3—stable subordinators. Section 3 introduces the definition of the generalized 
RL and Caputo type operators, and respectively. In Section 4 

we study important properties of the underlying stochastic processes associated 
with the generalized fractional operators. Then, the probabilistic solution to the 
equation o is addressed in Section 5. This section concludes with the application 
of these results to deduce the known solutions to fractional linear equations with 
the classical Caputo derivatives for (3 G (0,1). The last section presents the 

well-posedness to the mixed equations ([3]) and ©• 

2. Preliminaries 

2.1. Notation. Let N, C, and be the set of positive integers, the complex space 
and the d-dimensional Euclidean space, respectively. For any interval A, the stan¬ 
dard notation B{A), C'(H), Cb(A) and C'^(H) will denote the set of bounded Borel 
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measurable functions, continuous functions, bounded continuous functions and con¬ 
tinuously differentiable functions defined on respectively. If A is a closed set, 
C^{A) means the space of continuous functions differentiable up to the boundary. 
Notation || • || stands for the sup-norm ||/|| = sup^jg^ |/(x)| for / S B{A). If 
A = [a,b], notation Ca[a,b] and C\[a,b] denote the space of continuos functions 
vanishing at a and the space Ca[a,b] nC^[o, 6 ], respectively. 

As usual, Ep will refer to the Mittag-Leffter function of order ,5 > 0 defined by 




oo 


E 


r(j;9 + i)’ 


zee. 


whilst means the two-parameter Mittag-Leffter function: 


^Pi,132 (^) ^ ^ 


^.^r(j/3i + /32)’ 


z e C, /3i ,/32 > 0. 


For properties of these functions see, e.g., [7], [57]-[55]. 

Letters t and r are mainly used as space variables, and the letter s is used 
as a time variable. Bold letters, e.g. t and a, shall denote elements in for 
d > 2, and bold capital letters will stand for valued stochastic processes, e.g. 
T = {T(s) : s > 0}. Letters P and E are reserved for the probability and 
the mathematical expectation, respectively. Notation IF; 3 ((t, 7 ) means a /3—stable 
random variable (r.v.) with scaling parameter cr, skewness parameter 7 and location 
parameter zero. Its density function shall be denoted by Finally, for a 

given Feller semigroup on Cb{S), its resolvent operator at A > 0 is denoted 

by R\ and is defined as the Bochner integral (see e.g. m) 

poo 

Rxg := e-^^Ssgds, g€Cb{S). (10) 

Jo 

By taking A = 0 in (1101) . one obtains the potential operator which is denoted by 
Rog (whenever it exists). Additional superscripts shall be used to identify different 
resolvents and potential operators. 


2.2. Fractional diflferential operators. This section provides a quick summary 
of basic results concerning the classical Riemann-Liouville and Caputo fractional 
operators. For a detailed treatment refer, e.g., to [32]-I33] and references therein. 

Classical fractional differential operators are defined in terms of both the stan¬ 
dard differential operator (hereafter denoted by U™, m G N), and the Riemann- 
Liouville integral fractional operator /“_g for any a > 0, a G RU{—c»}, defined 

by 

1 /■* 

Ia+h(t) := / (t — s)°'~^h{s)ds, a > 0, t > a. 

r(a) 

For convention, /q_|_ refers to the identity operator. 

The left-sided Riemann-Liouville (RL) operator of order /3 > 0 (shortly the RL 
derivative) is defined as the left-inverse of the corresponding integral operator ([7], 
m)- To be precise, if /3 G R'*' and m = [/3] ( [•] denoting the ceiling function), 
then the RL derivative D^_^_ is defined by 

Df+/i(t) := D^LT+-‘'h{t), /3 > 0, /3 ^ N, t > a. 


( 11 ) 
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A sufficient condition for Z3f_|_ to be well-defined is to assume that / € ^"*[ 0 , 00 ), 
i.e., its derivatives of order m — 1 are absolutely continuous (see, e.g., Cl)- 

An alternative fractional differential operator is the left-sided Caputo operator 
(shortly the Caputo derivative): 

Di+Mt) ■= la+^D^hif), /3 > 0, /3 ^ N, t>a, (12) 

for which h requires the absolute integrability of its derivatives of order m = \ff\. 
It can be proved (see e.g. H) that both operators are related by the equality 

D^,^Mt) = Di^[h-T^_,[h-a]], (13) 

where Tm-i[h]a] denotes the Taylor expansion of order m — 1, centered at a, for 
the function h. Hence, in general 

i^f+/i(t) := D^C+-^h{t) # I^-PD^h{t) =■. 

unless the function h{t) along with its first m — 1 derivatives vanish at a-|- (or as 
t —>■ —00 for the case a = — 00 ). 


Remark 2.1. The left-sided derivatives have a direct counterpart to the right-sided 
versions (see previous references for details). However, we work everywhere in this 
paper only with the left-sided operators and their generalizations. The right-sided 
version of these results is a straightforward modification. 


2.3. Special case: /3 S (0,1). To solve fractional differential equations, we will be 
mostly interested in fractional derivatives of order /3 G (0,1). In this case equations 
m and (HID become 


A+Mi) = (/ (' - , t>a, 

J it- r)-t^h'(r)dr, t > a, 


respectively. Further, for smooth enough functions h (e.g. h in the Schwartz space), 
one obtains 


Dt+h{t) 


1 f* hit — r) — h{t) h{t) 

^+r(l-/l)(t-a)/5’ 


t > a, /3 £ (0,1), 
(14) 


and 

D^a+Mt) 


1 f* °‘h{t — r)—h{t) h{t) — h{a) 

n^Jo '■+r(l-/3)(t-a)/5’ 


t> a, 13 G (0,1), 
(15) 


for details see, e.g., Appendix in 

Thus, for these values of /3 the relationship between the Caputo and the RL 
derivates in m translates to 

D^+Mt) = D^+[^ - = D^+Ht) - > t > a. 

For smooth bounded integrable functions or functions that vanish at t = a (or as 
t -G —00 for the case a = — 00 ), the previous equality implies that the Caputo 
derivative and the RL derivative coincide. Its common value for a = — 00 , from 
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now on denoted by /dt^, is sometimes called the generator form of the fractional 
derivative of order jd € (0,1), [29]: 


dip 


hit) := DP_^^hit) = 


1 hit — r) — hit) 

Ti^Jo 


(16) 


2.4. Stable subordinators. Hereafter, we will always assume the existence of a 
probability space (fl, G,F) such that all the stochastic processes of our interest are 
defined on it. Notation means the completed natural filtration generated by a 
process X = {X(s)}s>o, he. := aiX^ : 0 < r < s). 

A ,5-stable subordinator for (3 G (0,1), is a real-valued stable Levy process TP = 
{T^{s)}s>o started at 0 almost surely (a.s.) with independent increments, T^(s) — 
T^(r) for any 0 < r < s, having the same distribution as the r.v. WpHs — r)^^P, 1), 
a totally skewed positive /3-stable r.v. with scale parameter cr = (s — see e.g. 

[H, [34]. 

This process has nondecreasing sample paths a.s., and it has only positive jumps. 
Moreover, this is time-homogeneous with respect to its natural filtration. Further, 
since ,5-stable processes are self-similar with index l/,5, the process {c^/^T^(s)}s>o 
has the same distribution as the process {T^(cs)}s>o for any positive constant c. 
Consequently, the transition probabilities i3) := P[r^(s) G i3|r^(0) = t] for 
any E G ,B(K) (the Borel sets of K) satisfy 


Psit,E) = s ^IP I wpis ^/^(r-t);l,l)dr, 

J E 

where wpi-] 1,1) is the density of a standard /3—stable r.v. Wgil, 1). The density 
in this case equals 


i(;; 3 (a:; 1,1) = ijR J exp |—TOa; — exp ^—7^/3^ | 


du, 


where 5R(z) means the real part of z G C (see Theorem 2.2.1 in [37)1. 

The infinitesimal generator of a ,5-stable subordinator, denoted by AP, is the 
generator of a jump-type Markov process of the form 


poo 

APhit)= / ihit + r) — hit))i'i3idr), hGDp, 

Jo 


(17) 


with a domain S,g and with the jump intensity given by the Levy measure v sup¬ 
ported in K-|.: 

= - ri-^)ri+P ^^- 

The last equality holds due to the identity r(a;) = (a; — l)r(a; — 1). 

For (3 G (0,1), we will say that the process T~^P = {r+^(s) : s > O} is an 
inverted (3—stable subordinator if —T~^P is a /3—stable subordinator. Thus, T~^P is 
a Markov process with non increasing sample paths a.s. and with the generator 


poo 

A~^Phit)= / ihit — r) — hit)) vpidr). 

Jo 


Notice that the relation 


wgi-t;aA) = w/3(t; cr,-1), 
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implies that T+^(s)—r+^(r) has the same distribution as the r.v. Wp ((s — — 1,0). 

Hence, the transition probabilities := P[T+^(s) G iii|T+^(0) = t] are 

given by 

E) = s-i/^ [ - r); 1, l)dr, E G H(K). (19) 

J E 


3. Generalized fractional operators 


Let = {T^(s) : s > 0} be a /3—stable subordinator with the generator and 
let = {T+^(s) : s > 0} be the corresponding inverted /3—stable subordinator 
with the generator . Then the operator —d^/dt^ in (HU) coincides with 

i.e. 

J/3 

-^Ht) = Dt^^hit) = = A+^h{t). 

As was shown in [22], an analogous probabilistic interpretation of the Caputo op¬ 
erators (resp. for any a G K and /3 G (0,1), can be obtained by 

interrupting (resp. killing ) on an attempt to cross the boundary point t = a. 
Moreover, this interruption procedure naturally yields an extension of the Caputo 
and the RL fractional derivatives. 

Namely, let generator of a decreasing Feller process taking values 

on (—oo, &], 6 G M, given by 



{h{t — r) — h{r)) p(t, r)dr, 


t <b, 


( 20 ) 


with a function v{t^r) satisfying the condition: 

(HO) the function v{t, r) is continuous as a function of two variables and continuously 
differentiable in the first variable. Furthermore, 


sup 

t 


rv{t,r)dr < oo, sup 

t 


dt 


v{t,r) 


dr < oo. 


and 

lim sup / rv{t,r)dr =Q. 

^^0 t J\r\<S 

Then, by Theorem 4.1 in [22], the process interrupted (and forced to land exactly 
at t = a) on the first attempt to cross the barrier point t = a is a Feller process on 
[a, b] and has the generator 


nt — G pCCi 

—Ll[+*/i(i) = / {h{t — r) — h{t))v{t,r)dr + {h{a) — h{t)) / v{t,r)dr, 
Jo J t—a 


t G [a, b] 
( 21 ) 


with the invariant core C^[a,b] and a domain denoted by Note that, since 

the process is decreasing, the interruption procedure effectively means stopping the 
process at the boundary. 

Moreover, if the process is also killed at the barrier point t = a (meaning analyt¬ 
ically to set h{a) = 0) then the corresponding Feller sub-Markov process on (a, b] 
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has the generator 


nt — a p<X) 

■.= / {h{t — r) — h{t))v{t,r)dr — h{t) / v{t,r)dr, t € {a,b], 
Jo J t—a 


( 22 ) 


with the invariant core C^[a, b] and a domain 

Definition 3.1. Let be a function satisfying condition (HO). The operator in 
(I 2 T]) will be called the Dt,-operator defined by v and will be denoted by (the 

sign — is introduced to comply with the standard notation of fractional derivatives). 
Similarly, the operator in (1221) will be called the D-operator defined by v and will 
be denoted by — 


Remark 3.1. Since we are interested in solutions to equations on finite intervals 
[a, &], we have applied the results in [ 22 ] to the case of stochastic processes taking 
values on (— 00 , 6 ] for some b gM.. 

By the standard theory of Feller processes, the domain of the generators 
and coincides with the images of their corresponding resolvent operators, 

j^a+*{v) (for any A > 0), respectively. Namely, 


£> 


a+* 



Ug(t) := g{t); g G C[a,b]j , 


and 


D 


M _ 

a+ 



Ug{t) := g G C'a[a,6]|. 


Moreover, the images of the resolvent operators are independent of A (for details 
see, e.g., 0 ). The operators in (1211) and (1221) can be thought of as the generalization 
of the Caputo and the Riemann-Liouville operators of order fi G (0,1), respectively. 

Particular cases of the D^— operators include the multi-term fractional opera¬ 
tors: 

- D^alMt) ■■= - J2ujfit)D^,XMt), 

2=1 


with non-negative functions > 0 and 


i/(t,r) 


d 

2=1 


1 


Even more generally, they include the case 

/ OO 

uj{s,t)Dl^^f^h{t) g[ds), 

-OO 

with 

v{t,r) = - 


(23) 


ds 


^‘"^"’^h(-/3(s,t))ri+/5(M)’ 


satisfying condition [HQ). Particular cases of (1^ have been studied e.g. in [^[T5] . 
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3.1. Classical fractional derivatives. For any /3 S (0,1), the fractional operators 
and are obtained as particular cases of D^- and H-operators, respectively. 
Namely, on smooth enough functions /i, 


if v{t,r) 


1 

r(_/3)r.i+/3 ’ 


P G (0,1) then 




-Dt+h{t). 


Thus, —is the generator of a Feller process on [a, h] which is obtained by 
stopping an inverted /3-stable subordinator in the attempt on crossing the 
boundary point t = a. The RL derivative is the generator obtained by 

killing at the boundary, see [35]. 


Remark 3.2. The probabilistic interpretation of fractional derivatives of order 
/3 G (1,2) was also analyzed in but this case will not be needed here. 


3.2. Fractional derivatives of position-dependent order. For a given func¬ 
tion /3 : R —>■ (0,1), define 

- T{-my*K‘y 

Lemma 3.1. If /3 : R (0,1) is a continuously differentiable function with values 
in a compact subset of (0,1), then the function defined in (l24l) satisfies condition 
(HO). 

Proof. Follows by the smoothness of the function /3 in a compact set of (0,1). □ 

Lemma lO allows us to define H*—operators by using They are denoted by 

—D^a+* — ~ ^a+l can be seen as generators of inverted stable-like processes 
m with jump density (1241) which are stopped at the boundary point t = a. These 
operators will be referred as Caputo-type operators of position dependent order. 
Analogously, the RL type operators defined by (l24l) will be denoted by —= 

7-)/5(t) 

^a-l- ■ 

Remark 3.3. Note that in the previous case the ’order’ of the derivative depends 
only on t. Some other extensions can be made by taking the function v depending 
on external variables. This case (which we analyze in detail in a forthcoming paper) 
allows us to deal with operators of the form 

( - ) hit, x), t>a,^e R^ 

where — denotes the Caputo type derivative acting on the variable t and 

depending on the variable x as a parameter; and — Ax ^ denotes the generator of 
a Feller process acting on the variable x and depending on the variable t as a 
parameter. 


4. Properties of the underlying stochastic processes 


In this section we study some facts about the underlying stochastic processes 
generated by the operators —and —£^1+. These results shall be used to obtain 
the explicit solutions to the linear equations involving 13,— and 13—operators. 

For a given function v satisfying condition (HO) and for t G (a, b], the notation 




-D 


a+*5 


t: 




and — D 


a-t-, 
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means the following: : s > 0} is the decreasing Feller process 

(started at t) generated by —as given in ([ 20 ]); : s > 0 } 

stands for the Feller process generated by —with the invariant core C^[a,b], 
which arises by interrupting on an attempt to cross the boundary point a; 

and : s > 0} denotes the Feller sub-Markov process generated 

by —D^a+ with the invariant core Ca[a, 5], which arises by killing at the 

point a. 

For t G [a, 6 ], notation refers to the first time the process (or the 

process leaves (a, 6 ], i.e. 

:= inf |s > 0 : ^ (a, 6] | = inf |s > 0 : ^ (a, 5] |, 




and, of course, = 0. Note that Ta^‘''^ is a stopping time with respect to 
Further, is also the first exit time from (a, b] of the killed process 

Let E) := P T^^''\s) G E \ = t be the transition probabilities 

of from tto E & B{—oo, b] during the interval [0, sj. Note that if E) 

denotes the corresponding transition probabilities of then 


p“+*M(r,i:;) = 


pt''’'\r,E), EGB{a,b] 

^ pt''''\r,{-oo,a\), E = {a}. 

Since the process is stopped at a, 

p:+*H(t, [a,t]) =p^M(t,(a,t])+p:+*M(t,{a}) = 1, 


r G (a, b]. 


and ps~^*^'^\a, E) = 1 whenever E G B[a,b] and a G E. Furthermore, iipB~^^''\r, E) 
denotes the transition probabilities of then ps~^^''\r, E) = pt^'^\r,E) for 

all r G {a,t] and E G B{a,b]. Moreover, 

it, (a, t]) = it, (a, tj) = 1 - pt^''^ it, i-oo, a]) < 1. 

The previous implies that 

pt^''Hr,E)=p-+<''\r,E)=p:+^''\r,E), r G (a,t], 
on Borel sets E of (a, 6 ]. 

Sometimes we will use the following additional assumptions concerning the func¬ 
tion v and the transition probabilities of the underlying process 

(HI) There exist e > 0 and <5 > 0 such that vit, r) > 5 > 0 for all t and all |r| < e. 

(H2) The transition probabilities of the process are absolutely continuous 

with respect to the Lebesgue measure (the transition densities will be denoted 
hy pt^''\r,y)). 

(H3) The transition density function pt^''\r, y) is continuously differentiable in the 
variable s. 


Remark 4.1. There exists several criteria in terms of v that ensure the validity of 
(iL 2 ) and (iL3), see, e.g., [15] . 


Lemma 4.1. Suppose the conditions (HO)-(Hl) 
stopping time is finite a.s. and E 


hold for a function v. Then, the 
< -foo uniformly on < G (a, 6 ]. 
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Further, the point a is regular in expectation for both operators —and 


i.e. E 




0 as 14 , a. 


Proof. The result follows by comparing the process with a compound Poisson 

process with Levy kernel v{r) = (5l{[g_g]|(r) for which the result holds. □ 

Remark 4.2. It is worth stressing that the assumption (HO) is weaker to the one 
used in Theorem 4.1 in [22] . 

Observe now that, since 13 g, decreasing process, the equivalence between 

the events | ^ ® > a | for t G (a, b] and all s > 0 implies 

> sj = P > a = [ p‘^+*^''\t,r)dr = 1 - / pt^''\t,r)dr, 

^ J (a,tl —00 


yielding the following result. 


Proposition 4.1. Suppose the conditions (H0)-(H3) hold. Then, the probability 
law of , denoted by (ds), is absolutely continuous with respect to Lebesgue 
measure for t G (a, b] and its density is given by 


= §;J Pt^''Kt^r)dr = J p+(‘')(t,r)dr. 


(25) 


We shall also need the joint distribution of and Ta^’'^ for any s > 0. 

Notice that for any a < r < t, 


r“+*(")(s) > r, > C = P rr*'"'(s) > r, > 




a+*(r/') 




Moreover, ^ < s implies 


y,a+*(,.)(^) ^ > aj = P [t“+*^‘'^(s) > r 


whilst for s < ^, 


T“+*(^)(s) > r, Tr'-''>{0 >a = I {t,w) ( / p“:r^^(n;,y)d 2 /) dw 


-^a+*{v 


a+*{u) I 


a+*(i') / 


= J / pt-s 


Therefore, 


■■= ^P < r, r*» < ^] = ^P [r“+*(‘')(s) > r, ^ 

yields the next result. 

Proposition 4.2. Suppose the conditions (H0)-(H3) hold. Then, for any s > 0 
and t G (a, 6], the joint distribution of the pair ^T“''’*^‘"^(s), Ta , denoted by 
Ps^a\dr,d^), has the density given by 

(r, 0 ■ = ^{s<i}Pt^''^ (t, P't-] (r, y)dy, a<r< t. (26) 


92 
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Remark 4.3. Since the processes j,a+*(v) rp+iv) before the 

first exit time Ta^''\ the equation (1^ provides the joint density of 
and (r“'''^^^(s), Ta for any s > 0 and for s < ^, respectively. 

By definition of the generator of a Feller process (see, e.g. [H], [2T]), if 5““'’*^'^^ 
is the semigroup of the stopped process ^j^gn u S 'S‘'a+* if) only if, 


— = lim 

“+ siO 


rra+*(i^) 

Ds U — U 


where the limit is in the sense of the norm in C[a, b]. Analogously, if is the 

semigroup of the killed process then u € if, and only if. 


— dI,''} u = lim 
+ s4.0 




where the limit is in the sense of the norm in Ca[a, b]. 

Let us now introduce an operator which will play and important role to charac¬ 
terize the domain of the generators and — 


For A > 0, define 

r+('^) 


:= E 




, te{a,b], 


for any (non constant) function g G B[a,b], and 


M+i‘'h(t) := E 




.. — As 


ds 


t e [a, 6], 


when git) = l(t) (the constant function 1). Then 


M 


H^) 




1 -E 




t,(u) 


(27) 


(28) 


(29) 


implying 

Further, 


E 




t,{u) 


= l-AM+i‘^).l(t). 


=0^7+1“^)-1(0, te[a,&], 

for any constant function equals to c (we shall use it mainly for the constant g{a)). 
Note that the stochastic continuity of the process implies that M^''lg(-) is 

continuous on [a, &]. Moreover, 


|M+M5(f)| < llffll sup E 


t^[a,b] 


r^,{v) 


Lemma 4.2. Suppose that v satisfies the conditions (H0)-(H3). Then 


E 




Jq ^ {^s J tG{a,b]-, 


and for any g G B[a,b] 


pz — u. pco 

^a\^9{t)= 9{t-r) e~^‘'pt^''\t,t-r)dsdr, t G {a,b]. 

Jo Jo 


(30) 


( 31 ) 
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Proof. Equality (|3n)l follows directly by using the density function of the r.v. 
as given in (l25)l . To prove ([ST]), observe that Fubini’s theorem allows one to 


rewrite as 




e"^®E 






ds. 


Using (l26)l . i.e., the joint density of the process for 

s < ^, yields 


PC 

= / 

^0 


-As 


Pt POO 

/ / 'i-{i>s}g{r)ipl’j^\r,^)d^dr ds 

Ja Jo 

J g{r)pt^''\t,r) (^ / pt-h^^y)dT?j d^drds 


-As 


g{r)p+^''\t,r) 


A r 

dlJ-. 


p^^'^\r,y)dy \ dj dr ds 


pt POO 

/ g{r)pt^''Ht,r) / yJ'^^''\-f)d-i dr ds 

Ja Jo 

g(r)p+('")(t,r) dr ds, 


pOO pt 


where the last equality holds as is the density function of the r.v. Ta^''\ 

This implies the result by another interchange in the order of integration and by a 
change of variable. □ 

Remark 4.4. Equality (l30l) can be written as 

/“OO 

^ — \ / ^ —As 


E 


= A / e 

^0 


Pt^''\t,r)dr) ds, t € {a,b]. 


(32) 


which follows by integration by parts. 

Let us now define the following space of functions: 

:= I M u{t) = • l(t) + d; t G [a,b], c, d G uj. (33) 

Lemma 4.3. Let z/ be a function satisfying conditions (HO)-(Hl). If A > 0, then 

1 


®i+* = <! Uo 


Ug{t) =9{a)j •l(t)) +M^^^^^g{t), g G C[a,b] 


and 


-^a+ — 


[wg Wg{t) = M+^‘'^5(t), g G Ca[a, &]|. 

Further, if also satisfies (H2)-(H3), then equalities (1501) and (1511) give explicit 
expressions for ^1 — ■ 1(<)^ and g{t), respectively. 

Proof. Let us take any u G ®i+*- Since is the generator of a Feller process 

on C[a,b], Theorem 1.1 in [5] implies the existence of a function g G C[a,b] such 
that u = 

u{t) = E 


that u = g. By definition of the resolvent and by Fubini’s theorem 


-XSg^j.a+*C')^s))ds 


UO 


1 


= E 

1 

O 

+ 


-XSg^j.a+H^)^s))ds, 
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where E 




< +00 by Lemma ICT 

Since the process is stopped at time g = 

Ta^’^\ Moreover, before time the processes and 

fore, 


g{a) for all s > 
coincide. There- 


■ pCO 



/ e-^^ds 

+ E 

/ e-^^g{T+^'^\s))ds 



Jo 




] 


E 

/ e-^"ds 

} +E 

/ e-^^g{T+^''\s))ds 


Jo 

J 

Jo 


= g{a)j (l - AM+M . l(t)) + M+^g{t), (34) 


as required. The characterization of the domain is similar to the previous case. 
Take any w € then there exists a function g G Ca[a, b] such that w = g. 

Hence, a similar procedure yields (IMll and g{a) = 0 implies 


^a+(v) 




9{t)- 


Finally, observe that under assumptions (H2)-(H3), Lemma [4.21 holds. □ Let 

us now see how the resolvents (and hence the domains) of the stopped and killed 
processes are related to. 


Lemma 4.4. Let be a function satisfying condition (HO). Suppose A > 0 and 
g G C[a,b]. Define g{t) = g{t) — g{a), then 


R' 




'>git) = • l(t), 


and 

i?y(^)g(t) = M^i^^git) - gia)M^i^^ • l(t). (35) 

In particular, R'^^''^g{t) belongs to both domains 1)^a+* 2)^+. 

Proof. Follows directly from the linearity of the operators and and 

by using that g{a) =0. □ 


Remark 4.5. Let us stress that Lemma [4.41 implies that coincides with 

the resolvent g only when the function g(a) = 0. Hence, only in this case 

belongs to the domain of both generators —and -Diy 

It is worth noting that Theorem I.l in [8] also guarantees that for g G C\a,b] 
and A > 0, the function Ug := g is the unique solution in the domain of 

-Dit to 

- D^yu(t) = Xu{t) - g{t), t G [a, b]. (36) 

Similarly, if 5 G Co [a, b], then w = g is the unique solution in the domain of 

-Di^i to 

- = Xw{t) - g{t), t G [a, b]. 
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5. Equations involving D— and D*- operators 

The probabilistic representation of solutions (in the generalized sense) to linear 
equations involving D- and -D*-operators will be studied in this section. 

5.1. Linear equations involving H-operators. Consider the problem of finding 
a continuous function u on [a, b] satisfying 

- D^^lw{t) = \w{t) - g{t), t€[a,b], w{a) = Wa, (37) 

for A > 0, g S B[a,b] and Wa = 0. Hereafter, we shall refer to (l37l) as the RL type 
problem {—Dl^^,X,g,Wa), wherein we shall always take Wa = 0. Similar notation 
will be used for the linear equation with the corresponding Caputo type operator: 
{-D[‘2*,\g,Wa) for any Wa e R. 

Definition 5.1. Let g G B\a, b] and A > 0. A function w G Ca[a, b] is said to solve 
the RL type problem X, g,0) as 

(i) a solution in the domain of the generator if w satisfies (1371) and belongs to the 
domain of the generator 

(ii) a generalized solution if for all sequence of functions gn G Ca[(i,b] such that 
sup„ ll^nll < oo uniformly on n, and liuin^oo dn —t 9 a-e., it holds that w{t) = 
lim„_).oo Wnit) for all t G [a, b], where Wn is the solution (in the domain of the 
generator) to the RL problem {—Dllfl, X, gn,0)- 

Remark 5.1. From this definition it follows that if there exists a generalized 
solution, then this is unique. 

Definition 5.2. The RL type equation (IWl) is well-posed in the generalized sense 
if it has a unique generalized solution. 

Well-posedness result for the RL type linear equation. 

Theorem 5.1. (Case A > 0) Let u be a function satisfying conditions (HO)-(Hl) 
and assume A > 0. 

(i) li g G Ca[a,b], then the linear problem (—A, 5,0) has a unique solution 

in the domain of the generator given hy w = g (the resolvent operator 

at A). 

(ii) For any g G B[a,b], the linear equation X, g,0) is well-posed in the 

generalized sense and the solution admits the stochastic representation 




- 

w{t) = E 

/ “ e-^’^g(T;^''\s)d 

5 ) 


Jo ^ 



Moreover, if additionally 1 / satisfies conditions (H2)-(H3), then 

pt — a poo 

wit) = / g{t — r) / e~^’^pf^''\t,t — r)dsdr. (39) 

Jo Jo 

(iii) If g G < 7 ( 0 , 6], then the solution to (l37l) belongs to direct 

sum of the domain of the generator and the space defined in (1551) 1. 
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Proof, (i) Take g G Ca[a,b]. Since g{a) = 0, w{a) = 0 and A > 0, Theorem 1.1 in 
[5] implies directly that w(t) = g{t) is the unique solution to (1571) belonging 

to the domain of the generator. Moreover, Lemma 14.31 implies 


n;(t)=M+Mg(t) = E 



(ii) Let us now take any function g G B[a,b]. As g does not necessarily belong to 
Ca[a, b], the resolvent operator no longer provides a solution to (1571) . However, using 
Definition IS.ll we will see that there exists a unique generalized solution. To do this, 
take any sequence gn G Ca[a,b] such that lim„_>oo = g a.e. and sup„ ||5n|| < oo 
uniformly on n. The procedure consists in finding the generalized solution as a 
limit of solutions to the equations 

- D^^lwn{t) = \wn{t) - gn{t), t G {a,b\, Wn{a) = 0. 


Since each gn G Ca[a,b], the previous case guarantees the existence of a unique 
solution Wn GD^^^I given by 


Wn{t) = E 



Using that ||gn|| is uniformly bounded, the dominated convergence theorem (DOT) 
implies 


lim Wn{t) = E 

n—>-oo 



=: w{t). 


Observe that the continuity of w{t) on (a, b] is a consequence of the continuity of 


the mapping < —>■ E 



The continuity at t = a holds by the regularity 


in expectation of (a consequence of assumption (HI)). Therefore, w G Ca[a, b] 
is a generalized solution to the linear equation (I37|) . Finally, the representation in 
(1551) follows directly from Lemma 1T51 


(iii) To prove that w G 0^+ ® 
in Lemma 14.41 to obtain 


a, A 


M. 


+(^ 


'9it)=K 


whenever g G C[a,b], we use the equality (|55l) 


where g(t) = git) — g(a), implying the result. 


□ 


Theorem 5.2. (Case A = 0) Theorem 15.11 with A = 0 is valid for the equation 

- D^alw{t) = -g{t), t G (a, &]; w{a) = 0. (40) 


Proof. It follows from the same arguments used for A > 0, so that we skip the 
details. Let us just notice that, since 




r \giT,^+^’^\s))\ds 

< Ibll sup E 

‘ a 

Jo 




Lemma l4 .1 1 implies that the potential operator corresponding to is bounded. 

Thus, Theorem 1.1’ in [8] ensures that the proof in 15.11 remains true for A > 0 if 
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one replaces the resolvent operator by the potential operator Also 

observe that tc G ®a+ ® whenever g G C[a,b] as w can be written as 


w{t) = + 5 '(a)E 

where g{t) := g{t) — g{a), for all t G [a, b]. 




ds 


□ 


5.2. Linear equations involving D*-operators. Let a G K and A > 0. Consider 
the problem of finding a function u G C[a, 6] solving 

- = Xu{t) - g{t), tG{a,b], u{a) = Ua, (41) 

for a given function g on [a, 6]. 

If u belongs to the domain of the generator then the linear equation 

(SD can be written in terms of the RL type operator as follows. Define 

w(t) := u{t) — Ua for all t G [a, b], then as = 0. 

Setting g{t) := g{t) — Xua, it follows that 

- = Xw{t) - g{t), tG(a, 6] w{a)=0. (42) 

Hence, u{t) = w{t)+Ua is a solution to the original problem if, and only if, w solves 
(1421) . The previous leads to the following definition. 

Definition 5.3. Let g G H[a, b] and A > 0. A function u in C[a, b] is said to solve 
the Caputo type problem X, g,Ua) as 

(i) a solution in the domain of the generator if w satisfies (j41ll and belongs to the 
domain of the generator 

(ii) a generalized solution if u(t) = Ua + w(t) for all t G [a,b], where w is the 
(possibly generalized) solution to the RL type problem {—D^J,X, g — Xua, 0). 

Remark 5.2. Definition 15.31 {ii), is given in terms of the RL type solution, but 
it can also be written in terms of the approximation of solutions belonging to the 
domain of the corresponding generator. 

Definition 5.4. The Caputo type equation m is well-posed in the generalized 
sense if it has a unique generalized solution and it depends continuously on the 
initial condition. 


Well-posedness result for the Caputo type linear equation. 

Theorem 5.3. (Case A > 0) Let i/ he a function satisfying conditions (HO)-(Hl) 
and suppose A > 0. 

(i) If g G C[a,b] and g{a) = Xua, then the linear equation {—D^J^, X, g,Ua) has 

a unique solution in the domain of the generator given by u = g (the 

resolvent operator at A). 

(ii) For any g G B[a,b] and Ua G K., the linear equation X, g,Ua) is well- 

posed in the generalized sense and the solution admits the stochastic repre¬ 
sentation 






u{t) = UaE 


+ E 

/ e-^‘^g(T+^'^\s)ds) 




JO ^ ^ 
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Moreover, if additionally u satisfies conditions (H2)-(H3), then 

pQO pt—a poo 

u{t) = Ua / + / g{t — r) / — r) ds dr, (44) 

Jo Jo Jo 

where pa^'^\s) denotes the density function of the r.v. Ta^''\ 

(hi) If g G C[a,b], then the solution to (HTl) belongs to 2)^^ © ^^‘'1 direct 

sum of the domain of the generator and the space defined in (l33l) l. 

Proof, (i) Since is the generator of a Feller process on C[a, b], for any g G 

C[a, b] the function u{t) := g{t) is the unique solution to (l36l) which belongs 

to the domain of the generator. This implies that u{a) = g{a) = g(a)/X. 

Hence, condition g{a) = Xua ensures that u satisfies the boundary condition in 
dm, as required. 


(ii) By Definition 15.31 u is the generalized solution if u{t) = w{t) + Ua, where w is 
the solution to the RL type problem (—, A, g(t) — Xua, 0) whose well-posedness 
is guaranteed by Theorem 15.II 

Moreover, the equality (1551) and the linearity of yield 


= Kx'^git) - Xu,M+^ -lit) = M+^g{t) - u, ( 


1 -E 


- — At, 


t,(u) 


where the last equality holds due to equation (l29l) . Thus, (l43|) is obtained by 
plugging the previous expression for w(t) into u{t) = w{t)+Ua- This representation 
implies directly the continuity on the initial condition Ua, as required for the well- 
posedness. Finally, the explicit solution (H51) follows from Lemma 14.21 


(iii) Assume now that g G C[a,b]. Since u{t) = g(t) — XuaM^'l ■ l(t) -|- Ua, 

by linearity one can rewrite it as 

- ff («) + 5 ( a )] (0 - • 1(0 + Ua 

= - 5(a)](0 + [g{a) - Xua]M^[''^ ■ l{t) + Ua- 

We then conclude that u G 

RT*‘'''\g - gia)] G 2)iy, 


and 

[g{a) + Xua]M+},^^-l{t)+Ua G 

□ 


Theorem 5.4. (Case A = 0) All assertions in Theorem 15.31 with A = 0 hold for 
the equation 

- A>i^*u(0 = -g(t), t G (a, 6]; u{a) = Ua- (45) 

Proof. Since the problem (1551) rewrites as 

- D^^lw{t) = g{t), t G (a, 6]; n;(a) = 0, (46) 

with w{t) = u(t) — Ua and g(t) = g{t), Theorem 15.21 gives the potential operator 
as the solution to (H51) for any g G Ca[a, 6]. Hence, the unique generalized 
solution to (1551) is given by u{t) = Ua + lim„_>oo gn(t) for any sequence gn 
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as in Definition 15.11 Consequently, the same arguments used for A > 0 remain 
valid. □ 


5.3. Fractional linear equations involving Caputo derivatives. Since Ca- 
puto derivatives are particular cases of the D*—operators, the solution to fractional 
linear equations with the Caputo derivative for fj G (0,1), is obtained by a 

direct application of the previous results. Namely, Theorem 15.31 implies that the 
problem 

D^+^u{t) = —Xu{t) + g{t), t e (a, b], u{a) = Ua G M, 
for a given g G B[a,b] and A > 0, has a unique generalized solution given by 





r 1 

r 

u{t) = UaE 


-f-E 

/ e-^^giT+^is))ds 




Jo 


where is the inverted /3—stable subordinator started at t. Moreover, substi¬ 
tuting (fT^ into formula (Hll) yields 


u{t) = Ua - a) J e ^ ^Wf 3 (^{t-a)s 1,1^^ ds-|- 

+ / git — r) / exp {—Asr^} 1, l)ds j dr. 

(48) 


Further, if (/(a) = Xua and g G C[a, b], then u belongs to the domain of the generator 


-D 


(C 

a+* * 


Corollary 5.1. Let t G (a, &] and A > 0. Then the Laplace transform of the first 
exit time from (a, b] for the inverted /3—stable subordinator started at t is given by 

E[e-^<'^] = Ep{-Xit-af), 

with Ep denoting the Mittag-Leffler function (see Preliminaries). Hence, 

1 1 / \ 

Epi —X{t — a)^) =—{t — a) J expi — Xs)s~'^~^wpyit — a)s~^^^',l,ljds. 

Proof. Follows as a consequence of formulas and (l4^ by taking g = 0 and 
Ua = 1 and by uniqueness of solutions. □ 

Remark 5.3. Alternatively, Corollary 15.II can also be obtained by using the iden¬ 
tity (see [37], Theorem 2.10.2) 


pOO 

PEpi-u)= / exp(-Mj/)?/“^“^/'^w;/3(j/“^'^'^; 1, l)dy. 

do 


(49) 


Moreover, using this identity we can see that 
solution given in ([S]). 


coincides with the well known 


6. Mixed linear equations 

In this section we study linear equations involving both the RL type and the 
Caputo type operators. The general setting will be explained first in and then 
we shall restrict ourselves to the simplest 2-dimensional case. This is done to avoid 
cumbersome calculations which nevertheless can be extended straightforward from 
the simple case analyzed here. 
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Let a = (ai,..., ad), b = {bi,... ,bd) G such that a < b. The space is 
assumed to be equipped with its natural partial order, the Pareto order, i.e. a < b 
means < bi for all i Notation [a, b] denotes the cartesian product 

[oi, 6i] X • • • X [ad, bd\ and t G [a, b] means ti G [a^, bi] for alH = 1,..., d. Let us 
denote by B[si, b] and C[aL, b] the space of bounded Borel measurable functions and 
continuous functions on [a, b], respectively, and by [a, b] the space of continuous 
functions on [a, b] with continuous first order partial derivatives on [a, b]. Similar 
notation will be used for (a, b] and (—oo,b]. 

Notation r“‘ means a vector r = (ri,... ,rd) G [a, b] having as its ith- 

coordinate. Since all the processes considered here have decreasing sample paths, 
we will be only interested in the boundary of (a, b] given by r“‘ for all z = 1 ,..., d. 
This subset will be denoted by 

d 

5a(a,b] := (J {r G [a,b] : r = }, 

i=l 

and the space of continuous functions on [a, b] vanishing at the boundary i9a(a, b] 
will be denoted by Ca[a, b]. 

For i = 1,..., d, let Vi be a function satisfying conditions (HO)-(Hl) and let 
t G (a, b]. The operator —tiD^ai+ represents the D-operator defined by Vi acting 
on the variable ti. For notational convenience set i' = [vi ,..., Vd) and define the 
mixed D—operator associated with the vector v as 


- ■L'a+ — - 2^ + > 


(50) 


Hence, the operator — is a sum of RL type operators each one acting on a 
different variable. Analogously, we define the mixed operators — D and — 
by using and respectively. 

Consider the RL type linear equation 


= Xw(t)-g{t), tG(a,b], 

w(t) = 0, tG9a(a, b]. 


(51) 


for a given function g € B[a., b] and A > 0. 

The operator — can be thought of as the generator of a Feller process on 
(a, b] obtained by killing the process generated by — on an attempt to cross 

the boundary 9a(a, b]. The killed process (started at t) will be denoted by = 

{t^(")(s):s>0}. 

Due to the independence, the Lie-Trotter theorem [111 |2T] implies that 
T^(‘^)(s) = (r“/+(‘^^)(s),... ,r“;+(‘^^)(s)), 

wherein each coordinate ig an independent process generated by —tiD^ai+- 

Let denote the semigroup of the process and let be the domain 

of its generator. Then, u G 2)2+ if) and only if, the limit 


s—^■O 


S 
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exists in the norm of Ca[a, b]. 

To solve m, let us introduce some definitions which will extend the ones used 
in the one-dimensional case. 


Definition 6.1. Let g € i3[a,b], and A > 0. A function w G C'a[a,b] is said to 
solve the RL type problem (— A, 5 , 0 ) as 

(i) a solution in the domain of the generator if w satisfies (1511) and belongs to 

(ii) a generalized solution if for all sequence of functions gn G C'a[a, b] such that 
sup„|| 5 „|| < oo uniformly on n, and gn ^ g a.e., it holds that w{t) = 
lim„_>oo iCn(t) for all t G [a, b], where Wn is the solution (in the domain of 
the generator) to the RL type problem (— A, 0). 

Remark 6.1. By definition, if there exists a generalized solution, then this is 
unique. 


For the sake of transparency, hereafter we restrict ourselves to the analysis for 
d = 2 and a = 0. Namely, let t = (< 1 ,^ 2 ), b = (&i,& 2 ) G t G [0,b]. Let us 
consider the equation 

- tiD^+^w(fl,f2) - t2D^+^'U'(fl,f2) = A'u(ti,<2) - 5(fi,f2), 

w(0, ^ 2 ) =0) = 0, 

where ti G (0, bi] for i = 1,2. 

Let pt^'''\ti,r) (resp. p^'^^'''\ti,r)) denote the transition density function of the 
process (resp. during the interval [0, s]. If is the first exit 

time from (0,&i] of the process (started at U), then the first exit time from 

(0, b] = (0, bi ] X (0, 62 ] of the process denoted by equals 

To ' ' = mm I Tp ’ : 1 = 1 , 2 ) . 

Due to the independence between the coordinates of the process , its transition 
density function, denoted by (t, r), satisfies 

i=l 

yielding the following result. 


Lemma 6.1. Let t = (fi,< 2 ) G (0,&i] x ( 0 ,& 2 ]- Suppose (HO)-(Hl) hold for both 
functions vi and V 2 - Then, 

(i) The boundary points ( 6 ,^ 2 ) € for all t 2 G [ 0 , 62 ), and (ti,0) G for all 
ti G [ 0 , 61 ), are regular in expectation for both operators —Dg*)!^ and — 


Moreover, E 




< -boo uniformly on t. 

(ii) If additionally each izi satisfies assumptions (H2)-(H3) and pQ^'^\ds) denotes 
the probability law of Tq’^'^\ then its density function fJ-o^'^\s) is given by 




s) = 


Mo 


X'^i) 


is) 


n+M 


(is 




(s) [ pt''''"Hti,r), s > 0 . 
Jo 
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(iii) Further, assuming again that each Vi also satisfies (H2)-(H3), the joint distri¬ 
bution of the pair denoted by (dr, the density 

f p'^^7\n,y)dy + 

Jo 

+ Js]o‘'^\n,Opt^''^Ht2,r2) [ P^i7^(^2,y)dy, 

do 

for all 0 < s < and r = (ri,r 2 ) G (0,t]. 

Proof, (j) The regularity in expectation of the boundary 9o(0, b] is a consequence 
of assumption (HI) and the Lyapunov method (same reasoning as in [22) 1 applied 
to the Lyapunov function 

da;(ti,t2) = ^lH2^ wi,a;2G(0,1). 

The finite expectation of is a consequence of the finite expectation of each 

T 

(ii) This is a generalization of Proposition id.ll and follows directly by differentiating 


P 


>s =P P To 


t 2 ,(i^ 2 ) 


> s 


with respect to s. Notice the use of the independence assumption in the previous 
equality. 

(iii) This is a generalization of Proposition and is obtained by differentiating 

2 


p 0 -|-(i 2 ) , 


i=l 


is) > > e = HP («) > > 0 


with respect to ri,r 2 and □ 

Let us now generalize the definitions given in (l?7l) and (1^51) . For A > 0 and 
g G i?[0, b] define 


M+(;)g(t) := E 


t,lu) 


e-^^giTt^''\s))ds 


, tG(0,b], 


and 


M+Ml(t) := E 


t,(.) 


e ^^ds 


, tG[0,b]. 


Note that Mq ^_('^ 5 (-) is continuous on [0, b] and |Mq ^_('^ 5 (t)| < HsH supt^jQ ^] E 
Moreover, 




• l(t) = i (l 


-E 




t.c,-) 




implying 


E 


,-At, 




= 1 - AM, 


0,A 


l(t), 


and yielding the next generalization of Lemma 


Lemma 6.2. Let t = (^ 1 ,^ 2 ) G (0, b] and A > 0. Suppose that satisfies condi¬ 
tions (HO)-(Hl), for i = 1,2. Then, 


E 






= E 


e-^^0 1 




l)<^t 2 ,(.' 2 )| 


E 


e “ l|^t2,(-2)<.rti,(-i)| 
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If additionally Vi satisfies (H2)-(H3) for each j = 1, 2, then 

(■*2 


E 


,-At, 




I 


^ —As 




Further, 

Mq A'^5(t) = / [ 9{ti-ri,t2-r2) [ {t 2 ,t 2 -r 2 ) ds dr 2 dri. 

Jo Jo Jo 

(52) 

Proof. Similar to the proof of Lemma 14.21 but using the density function of the 
r.v. and the joint distribution of the pair given in Lemma 

16.11 □ 
Well-posedness result for the RL type linear equation. 


Theorem 6.1. (Case A > 0) Let v = {i'i,i' 2 ) be a vector such that each i/i is 
a function satisfying conditions (HO)-(Hl). Suppose A > 0 and (ti,t 2 ) G (0,6i] x 
(0,62]. 

(i) If 5 G C'o[0,b], then the equation (— A, 5 , 0) has a unique 

solution in the domain of the generator given by w = the resolvent 

operator of the process 

(ii) For any g G R[0,b], the mixed equation (— tif5o+^ ~ t 2 -C>o+\A, 5 ,0) is well- 
posed in the generalized sense and the solution admits the stochastic repre¬ 
sentation 


w{ti,t2) = E 


^(tl42),(«') 


-Xsg ^ 


ds 


(53) 


Moreover, if additionally each i/i, i = 1,2, satisfies conditions (H2)-(H3), then 
w{ti,t 2 ) takes the explicit form in (15^ . 


Proof, (i) Follows from Theorem 1.1 in [5] as in the one-dimensional case. 

(ii) If g G R[0, b], the solution is obtained as a limit of solutions R°^^^^(/„(t) 
in the domain of the generator —where the sequence of functions {gn}n>i 
satisfies the conditions of Definition 16.11 Finally, Lemma 16.21 provides the explicit 
representation of the solution w in terms of transition densities. □ 


Theorem 6.2. (Case A = 0) All assertions in Theorem 16.II are valid for A = 0. 

Proof. The arguments in the proof of Theorem 16 .1 1 remain valid for the case A = 0 
replacing the resolvent by the corresponding potential operator □ 

Finally, we shall analyze the mixed linear equation which involves both the RL type 
and the Caputo type operator: 

- tiD^o+^u{ti,t2) - t2Do+lu{ti,t2) = Xu{ti,t2) - g{ti,t2), (ti,t2) G (0,61] X (0,62], 

U{0,t2) = 0, t 2 G[ 0 , 62 ] 

u(ti,0)= tiG( 0 , 61 ], 


( 54 ) 
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for a given function (j) € Co[0,6i]. This equation will be referred to as the mixed 
linear problem (- \ 9, (f) ■ 

Denote by := Feller process (with values on 

(0, 6 i] X [ 0 ,& 2 ]) generated by the operator — tiD^+^ ~ This process is 

obtained from a process := by either killing it whether 

the first coordinate attempt to cross the boundary point ti = 0 , or by stopping it 
if the second coordinate does the same with the boundary point t 2 = 0. As before, 
denotes the first exit time from ( 0 ,&i] x ( 0 , & 2 ]- 
In order to solve the mixed equation (1541) . we shall rewrite it as a linear equation 
involving only RL type operators. Namely, let ip £ x [ 0 ,^ 2 ]) be a function 

satisfying the boundary conditions in (|54ll . Define w{t) := u(t) — ip{t) for any 
t = (^ 1 ,^ 2 ) S [0,b]. Observe that, by definition, w vanishes at the boundary 
5o[0,b]. 

If u and Ip belong to the domain of the generator — — t 2 -Do+i, then 


(- 


ti^0+ 



tl^0 + 


t 2 ^ 


02 )\ 
0+* J 



ti^0+ 


+ 


* 2 ^ 


o+* 



= Xu — 



ti^0+ + 


* 2 ^ 


o+* 



Xw-g, 


with g := g — Xip — tiDQ'^'^ip — Due to the properties satisfied by ip, the 

function w satisfies that — ( 2 ^ 0 +*'^ = “ * 2 ^ 0 +^^ = 0 at the boundary 9o(0, b]- 
Consequently, the solution u to (1541) can be written as u = w + ip, where w is 
the solution to the corresponding RL type equation. This leads us to the next 
definition. 


Definition 6.2. Let g £ i?[0,b], A > 0, and (p £ Cop, 61 ]. A function u £ C[0,b] 
is said to solve the mixed linear problem (— A, 5 , (p) as 

(i) a solution in the domain of the generator if u satisfies (1541) and belongs to the 

domain of the generator — — t^Dp+i; 

(ii) a generalized solution if for any function ip in the domain of — — * 2 - 00 +] 

such that ip{0, •) = 0 and ip{-, 0 ) = <p{-), then u = uj + ip, where w is a solution 
(possibly generalized) to the RL type problem 

(- - rf;\A,5,o), 

with g := g - Xip - “ t 2 -Do+i' 0 - 

Remark 6.2. By definition, it seems that a generalized solution depends on the 
function ip, the next result we will show that this is actually independent of ip. 

Remark 6.3. The concept of a generalized solution can also be given as a limit of 
solutions in the domain of the generator. 

Well-posedness result for the mixed linear equation. 

Theorem 6.3. (Case A > 0) Let v = {i/i, V 2 ) be a vector such that each Vi is a 
function satisfying conditions (HO)-(Hl). Suppose A > 0 and (p £ Co[0, 61 ]. 

(i) If 5 G C[0,b] satisfies g(0,-) = 0 and g(-,0) = Xfp(-), then the mixed equa¬ 
tion (— — t 2 ^o+I, X, g, p) has a unique solution in the domain of 
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the generator given by it = the resolvent operator of the process 

^ 2 ^ 0 +(!^ i ) fj,0+*(v2.) ^ 

(ii) For any g G B[0,b], the mixed linear equation (— tiD^+ — t 2 -Do+i, A; ffi </^) 
is well-posed in the generalized sense and the solution admits the stochastic 
representation 


u{ti,t2) = E 


e ° 




rpO+{vi) 


^0 


l/_*2.(‘'2) 


,(- 1 ) 


} 


+ E 



(55) 


Moreover, if additionally each i = 1,2, satisfies conditions (H2)-(H3), then 
the solution rewrites as 

pii /‘oo 

u{ti,t2) = / - ri) ds dr-f 

Jo Jo 

nt\ nt2 poo 

+ / g{ti - ri,t 2 - r 2 ) e~^‘^pf^''^^ti,ti-ri)pf^''^'^{t 2 ,t 2 -r 2 )ds dr 2 dri. 

Jo Jo Jo 

(56) 

Proof, (i) As before, we apply Theorem 1.1 in [8] to the generator — — ta-Do+i • 

Therefore, if (/ is a continuous function on [0, bi] x [0, 62 ] such that g(0, •) = 0, 
then the function u(ti,t 2 ) = ^ 2 ) solves the equation (IMl) without any 

boundary condition. Further, a simple calculation shows that 

u(ti, 0) = 0) = g(ti, 0)/A, 

which implies that, under condition g{-,0) = A())(-), the function u solves the 
problem (IMl) . 

(ii) For the general case g G 5(0, b], take a function -0 satisfying the conditions 
of Definition 16.21 and set w := u — ip. Since w vanishes at the boundary 9o(0, b]. 
Theorem 16.11 yields 


w{t) = E 


,(tl,‘2),(*') 


^ —As 


~9 , T°+(‘'=)(s)) ds 


with g = g — Xip — { tiD^+ + i 2 'C*o+*) V'- Hence w{t) = I — II, where 

,(tl.‘2).(‘') 

e-^^g ( t °+("^)( s ) , ds 




I- II := E 


-E 


10 


„ — Xs 


(A+ tA7'’+ > TA’'^\s))ds 


Using that ip belongs to the domain of the generator — tiOA ~ * 2 ^ 0 +* > Doob’s 
stopping theorem applied to the martingale 
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(57) 


+ 




A + ti-Dp 


i’^i) 


+ 




with the stopping time 7 ^* 1 jjj^pijeg 


II = 1p{ti,t2) - E 




^ ^2i0+(l^l) j'^(tl,*2),(!^)^ 'Ji 0 + (!22) j'^(tl,t2),(!2)^^ 


which, in turn, yields ([55|l as u = w + ip and ip{-, 0) = ^(•). 

Finally, the second term in (I56p is a consequence of Lemma |6.2l whilst the first 
term is obtained by conditioning first on Tq^A'^^) then by using the joint density 


of the pair (s), 




□ 


Theorem 6.4. (Case A = 0) All the assertions in Theorem 16.31 are valid for the 
case A = 0. 


Proof. For functions g G Co[0,b], the arguments in the proof of Theorem 16.31 
remain valid using the potential operator IHq ^ ^ instead of the resolvent operator 
^o+(i/)*^ For the general case g G i3[0,b], the martingale (1571) should be replaced 
by the corresponding martingale with A = 0. □ 

Remark 6.4. As an application of Theorem 16.31 one obtains that for t = (ti,t 2 ), 
the function 

u{t) ——t2 j — r) J 1 ,ui/3 1 ,ds dr + 

+ J J g{ti - ri,t2 - r2) J Wa (r2S~^^°‘-,l,l^ dsdr2dri 

is the generalized solution to the mixed fractional linear equation 

- tiD^+u{ti,t2) - t2Do+^u{ti,t2) = \u(ti,t2) - g{tl,t2), (tl,t2) G (0,b], 

u(0,t2) = 0, i2G[0,62] 

u(ti, 0 ) = (j){ti) tiG(0,5i], 


for a given function (f G Co[0,5i] and P,a £ (0,1). Let us recall that — tiDg^ 
and — t 2 Do+* stand for the classical RL and Caputo derivatives of order /3 and a, 
respectively; and Wjs and Wa denote j5— and a—stable densities, respectively (see 
Preliminaries). 

The solution u belongs to the domain of the generator only when g G C[0,b], 
g(-, 0 ) = and 5 ( 0 ,-) = 0 . 
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